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■ Abstract 

OO 



We study theoretically and numerically the effects of the linear velocity 



■ field V = voyx on the irreversible reaction A + i? — > 0. Assuming homo- 

geneous initial conditions for the two species, with equal initial densities, 



we demonstrate the presence of a crossover time tc ~ ^- ^ ^ ^' 
the kinetics are unaffected by the shear and we retain both the effect of 



Q , species segregation (for d < 4) and the density decay rate At where 

O ■ 

a = min(^,l). We calculate the amplitude A to leading order in a small 
^ ' density expansion for 2 < d < 4, and give bounds in d = 4. However, for 

t ^ Vq , the critical dimension for anomalous kinetics is reduced to dc = 2, 
with the density decay rate Bt~^ holding for d > 2. Bounds are calculated 
for the amplitude B in d = 2, which depend on the velocity gradient vq 
and the (equal) diffusion constants D. We also briefly consider the case of 
a non-linear shear flow, where we give a more general form for the crossover 
time tc. Finally, we perform numerical simulations for a linear shear flow in 
d = 2 with results in agreement with theoretical predictions. 

PACS Numbers: 02.50. -r, 05.40. +j, 47.70. -n, 82.20. -w. 
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1 Introduction 



Recently there has been considerable activity in the field of diffusion limited chemi- 
cal reactions (see [1-9] and references therein). Many studies have concentrated on 
the effects of density fluctuations, especially in the one and two species reactions 
A + A^^ and A + i? ^ 0. In the time dependent case it is well known that for suf- 
flciently low spatial dimensions these fluctuations alter the kinetics. For example, 
in the two species reaction with equal initial concentrations, the densities decay 
asymptotically as for (i<4|]l], ^|^-a slower rate than the mean fleld 
result. This is closely related to the phenomenon of segregation, where the species 
separate into A and B rich zones at large times (for d < 4). These effects are 
essentially due to inadequate diffusive mixing of the reactants in low dimensions. 
This allows the initial density fluctuations to persist, leading to the formation of 
segregated zones. 

However, exposure of the system to a shear flow will lead to a modiflcation 
of these kinetics, as the velocity gradient will allow the chemicals to mix more 
efficiently. Hence, we might expect that the critical dimension for the onset of 
anomalous kinetics (and the appearance of segregated zones) might be lowered. 
In this paper we attempt to confirm these intuitive expectations by focusing on a 
two species reaction in a linear (Couette) shear fiow, in the case where both species 
have the same diffusion constant D, and the same initial density no. Previous work 
on similar problems includes the study of shear forces on a binary fluid mixture at 
criticality (see |]TD[ and references therein). More recently studies have been made 
of reaction-diffusion systems where particles with a same species exclusion rule were 
subject to a drift [jn], |r^, [TB[. This exclusion rule meant that particles of the same 
species were forbidden from occupying the same lattice site. New exponents for 
the asymptotic density decays were reported, with a theoretical justification based 
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on the Burgers equation. Note that our model differs from these cases, both by the 
nature of the velocity flow (shear not drift), and by the absence of any exclusion 
rules. 

Although very little work has been done on fluctuations in reaction-diffusion- 



shear systems, a considerable amount is known about diffusion in shear flows [IJ, |15 



T6| , p!7| . Exact solutions have been given for the positional probability distribution 



of a Brownian particle released in a linear velocity field. These results will be 
employed in the analysis of the following sections. We also mention a paper by Yee 
[|T^ , where a method was given to obtain an approximate analytic solution to a 
general reaction-convection-diffusion equation. However, this technique is simply a 
perturbative method for solving the mean field equations, and as such it takes no 
account of the microscopic density fluctuations. So instead we employ a different 
method based on a mapping of the microscopic dynamics onto a quantum field 
theory 0, |], [l^, This allows for a systematic treatment of the density 



fluctuations using diagrammatic perturbation theory. In this way we are able 
to go beyond the traditional rate equations approach, by using the calculational 
framework of either an effective or a full quantum field theory. 

In brief we find, for a linear shear flow with velocity gradient Vq, a crossover time 
of order Vq^ - at times t <^ Vq^ the behaviour of the system is essentially unaffected 
by the shear. However for t ^ Vq^, the critical dimension for the system is reduced 
from dc = 4: to dc = 2, and hence the mean field decay exponent holds in all 
physically realisable dimensions. In addition, we note in passing that the problem 
of the single species reaction A + A ^ in a shear flow is unlikely to be interesting, 
as its critical dimension even without shear is already (ic = 2 |1|, 0, |]. 

Finally, we give an outline of the layout of this paper. In section 2 we give 
simple arguments for the critical dimension and crossover times for these systems. 
The mapping to a quantum field theory is carried out in section 3, and in section 



4 we perform density calculations in the different regimes, using the field theory 
formalism. Numerical simulation results are given in section 5, and we present our 
conclusions in section 6. 



2 Decay Rates and Crossover Times 

We shall first consider some simple arguments which determine the density decay 
rates and crossover times induced by the presence of a shear flow. If we neglect 
the role played by fluctuations, then the mean field result for the densities gives an 
asymptotic decay. However by using a variant of an argument first put forward 
by Toussaint and Wilczek , we can understand how the density fluctuations alter 
this result. At t = 0, these fluctuations ensure that 

|iVa-iV6|t=0~K^')^ (1) 

where Na,Ni, are the number of A,B particles within a volume L'^. Consequently 
the initial density difference satisfies: 

(2) 

In the absence of a shear flow, after a time t, typically only the initial density 
excess will remain in a volume of size (Dt)^ , as all the other particles in that 
region will have mutually annihilated. Hence the number of particles remaining in 
volume (Dty^'^ at time t is {no{Dt)^)^ (using equation (|l|)) and we thus obtain a 
t~i density decay |^, |^. However, in our case, the presence of the shear flow 
means that length scales parallel to the velocity flow (in the x direction) increase at 
a different rate to length scales in perpendicular directions. If we now specialise to 
the case of linear shear, then it is known (from exact solutions for random walkers 



in linear velocity flows ||14| , p!5| , \W^ ) that this characteristic length scale grows as 



{Dt[l + i(fot)^])2. Consequently, after a time t, only the initial density excess will 



remain in a volume of size (Dt)''^^ {Dt)2[l + |('i^o^)^]^- As a result we expect the 
densities to decay asymptotically as 



1 + 



(3) 



Thus for t <^ fp"^, we retain the t~i decay (for d < 4), with the mean field 
exponent applying for > 4. However for times very much larger than the crossover 
time tc ^ Vq^ ^ we have a different regime: 



(d+2) 
t 4 t) 



for < 2 



a^h^ r (4) 
t-^ for d>2. 

Note that c? = 2 is the lowest possible dimension for the geometry of our system, 
so in practice the mean field decay exponent is always retained. 

We now briefly consider the case of non-linear shear - a considerably more 
complicated situation. However, we can give some simple arguments which reveal 
the crossover time where we expect the shear to begin to alter the dynamics. 
Consider, for example, the non-linear flow v = fo|y|"sgn(y)x studied in The 
shear flow will disrupt a segregated zone around y = yo when the typical zone 
length scale (Dtc)^ is comparable to the distance over which the top and bottom 
of the segregated zone are sheared apart. The time tc at which this happens is 
given by 

VoUiDQ"^ + yoT - votcy^ ~ {DtS^. (5) 



We can rewrite this relation for as 



votcyo 



1 + 



yo 



(6) 



For the case of a linear shear flow (n = 1), we have recovered our earlier result of 
a crossover time tc ~ Vq^- However for n 7^ 1 and (Dtc)^ ^ yo, we can expand the 
bracket in the above equation. This leads to a crossover time 



tc 



vonyo 



_d_ 

dyo 



1 -1 



(7) 



valid when 



D2 1 \2 D 
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3 The Field Theory Approach 

In order to perform more quantitative calculations for the densities, we need a 
systematic way of including the effect of microscopic density fluctuations. One 
way in which this can be achieved is by mapping the microscopic dynamics onto 
a quantum field theory. The first step in this process is to write down a master 
equation for the microscopic dynamics. This can then be recast in a second quan- 
tised formalism, which may in turn be mapped onto a path integral. All of these 
steps have been described in detail elsewhere ^ |, [T^, so we simply give 
the resulting field theoretic action. Defining the fields (f) and ip in terms of the 
continuous c- number fields a and 6, where (p = |(a + 6) and ^/^ = |(a — 6), we have: 

S = j dxdy<f-'^zdt \2i){dt + v{y)d,^ - DV^)iIj + 2(f){dt + v{y)d,^ - DV^)(j) 

+2A0(02 - + A(02 - (02 _ _ 2no05(t)] (9) 

where and (f) are the response fields. The vertices for this field theory are shown 
in figure 1, where the propagators for the field are solid lines, and the ip propa- 
gators are dotted lines. Notice that if we neglect the quartic terms in the action, 
and integrate out the response fields, then we recover the classical (mean field) 
equations: 

f.„(,)|^^.VV (10, 

^ + ^(^)|^ = ^^'^ - - + ^o^(^)- (11) 
To derive the form of the Green functions we again specialise to a linear velocity 



« 1. (8) 
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field, i.e. v{y) = Voy. The equation for tlie Green functions G = G^^ = G^^ is: 

2 + voy^ - DV'^ G = 6{x- x')6{y - y')6{z - z')6{t - t'). (12) 



Tfiis equation lias been solved elsewhere [T^, [15], [T6|, |T^, in the context of the 
position probability distribution for a random walker released from the origin in a 
linear shear flow: 



G^iA.D(t^f))-^[-^^-^^\ X (13) 



exp 



3(a: -x'- \v,{y + y'){t - t')f {y - y'f (z - z'f 



D{t- f) [{vo{t - t')Y + 12] AD{t - t') AD{t - t') J 

Note that for fo(t — t') ^ 1 the Green function is the essentially same as for a 
shearless system (with fo = 0). Physically this result corresponds to the dominance 
of diffusion over shear on short enough time scales. 

Finally, we can use the field theory to define two classical (tree level) quantities. 
The first of these is the classical density (0)c/, which is the sum of tree diagrams 
contributing to (0) (see figure 2, where (0)c/ is represented by a wavy solid line). 
This sum can easily be evaluated 0, |^, giving 

the same result as would be found by solving the mean field rate equations. The 
second classical quantity is the response function, which is defined to be the sum 
of all possible tree diagrams connected to a single propagator. The nature of the 
vertices ensures that the ip response function is just the propagator, whereas the 
response function is given by the diagrammatic sum shown in figure 3 (where it 
is represented by a thick solid line). 
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4 Field Theory Density Calculations 

Using the formahsm presented in section 3, we are now in a position to calculate the 
densities in the two distinct time regimes. We consider first the case of 2 < < 4 
and Vot -C 1, where an effective field theory can be developed. However such a 
theory turns out to be almost identical to that previously used by Lee and Cardy 
to study the same two species reaction, but without shear (see for details). The 
only fundamental difference lies in the slightly modified form of the Green functions 
(|T3|). However, they are still sufficiently similar to those of Lee and Cardy that 
our results, to leading order in a small density expansion, reproduce those of [0, 
and are independent of vq. Hence, quoting from 0, we have, for 2 < (i < 4 and 
(noA)"^ -C t < Vq^: 



leading order in a small density expansion. Note that we require t ^ (noA)^^ so 
that the coarse-graining required for the calculation of the initial term A is valid. 
Strictly for this and subsequent results, we must also have Dt ^ ng'^^'^ - i.e. the 
particles must have had time to "find" each other by diffusion. For c? = 4 we have 
the upper bound 



where Ae// is an effective coupling constant found by summing the diagrams in 
figure 5, giving for t>ot -C 1: 



Here A is a cutoff which is needed to keep the loop integrals finite. The lower 
bound for the densities in c? = 4 is given by 





(16) 




(17) 



{b)iowi 



-1 



(18) 
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Finally, it is important to note that these calculations are valid only in the regime 
(tiqA)^^ <^ t Vq^ and Dt ^ Uq '^^'^ . Hence, for small enough and large enough 
Vq, these conditions will not be satisfied and the above results will not be applicable. 

4.1 v^i <C 1 and d = 2 

Power counting on the full field theoretic action (P) gives c? = 2 as a critical dimen- 
sion for the system (when v^t ^1). In this case (of the lowest physically possible 
dimension), we must consider the full theory, as given by the action (^). The renor- 
malisation is similar to that previously developed in 0, where more details 
may be found. In particular the field theory remains simple in that diagrams cannot 
be drawn which dress the propagators. Consequently the bare propagators are the 
full propagators and both f o and D are not renormalised. Furthermore the simpler 
form of G for v^t <^ 1 again ensures that the results for the (0^) and ('?/^^) loops 
are unchanged from the shearless case 0. Thus the primitively divergent vertex 
function X{k,t) (figure 6), and resulting (3 function, remain unaltered. Note that 
we have now rescaled the couplings to absorb the diffusion constant D: X = X/D. 
The Callan-Symanzik equation for the densities is slightly modified to read 



- (9 d d d d 

21-^ - 2vo— - dno^ dA— + Pignj^— + d 

at ovq ouq oA ogR 



n{t,gR,no,vo,A) = 0, 

(19) 

where vq = vq/ D, i = Dt, and g^ is the dimensionless renormalised coupling. The 
solution can be found by the method of characteristics 

n{t, gR, no, Vq, A) = {K^t)-^n{K-'^, gR, no, #o, A), (20) 

where no = no(K^t)^, A = A(K^f)^, and vq = vo{K,'^t). Furthermore, at large 
enough times (but still such that Vot <^ 1), the running coupling gR goes to zero 
as (lnf)~^ for d = 2 The leading order result is given by (|15|) , and following 
[0, we make the assumptions that higher order terms are both independent of no, 
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Vq (and thus of no, "^0)5 and that they diverge no more quickly than As, for large 
A. Consequently, if these assumptions are valid, then the densities in = 2 will be 
given by expression (|15]), with corrections which are suppressed by at least a factor 
of (Int)-^ 

4.2 vot > 1 

At truly asymptotic times t ^ Vq^, we can use the Green functions G = G.^^ = G^^ 
(|T3D to motivate a new assignment of dimensions for the parameters appearing in 
the action (P). Firstly we can see from (|13D that [x\ ~ [Dv^t^]^ . If we ignore the 
d1 terms in the action (corresponding to the neglect of diffusive motion in the x 
direction) and assign [vq\ ~ , then we can give the following naive dimensions: 

[x] ~ k-^ [y], [z] ~ fc-i [t\ ~ A;-2 [A] ~ k-'' [0], [^] ~ [0], [^] ~ k''+\ 

(21) 

This suggests that a full field theory analysis using the action (^), with subsequent 
renormalisation, becomes necessary only at dimension c? = 0. 

Consequently, following , we must now construct a new effective field theory, 
valid for d > 0. The first step in this process is to determine which initial parame- 
ters are relevant. For an initial term of the type (A*^™''"ym!n!)0'"-?/'"|j=o, we must 
therefore consider the dimensions of the coupling [A^+^A*^™'")] ~ ^-d{m+n)+d+2 ^ 
The power of A follows from considering the number of vertices needed to attach the 
initial term to a given diagram. These terms will be relevant when d < 2/(m+n— 1). 
Hence, if we have m + n = 1, then such an initial term is relevant for all d. The 
case of m = 1 corresponds to the initial density, whereas the generation of an n = 1 
term is forbidden by the invariance of the system under a transformation exchang- 
ing A <^ B, i.e. (0,0,-0,'?/') — i> {(f),(j),—ilj,—ilj). The only other important initial 
terms are those with m + n = 2, which are marginal in ci = 2. In fact we need 
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only consider an extra initial term of the type (A*^°'^Y2)V^^, as A^^'^^ is forbidden 
by symmetry, and A^^'^^ is suppressed (it can only act as a source of noise through 
a response function, which we assume to be heavily damped, as in 0). So, for 
d < 2, we are led to the construction of an effective field theory with an extra initial 
term, whereas for d > 2 all initial terms (except the initial density) are irrelevant 
and hence the rate equation approach can be employed. In what follows we shall 
develop the effective field theory only in ci = 2, as the system cannot be realised in 
a lower dimension. 

Turning to the calculation of A*^"'^^ = A, we again need to sum the set of 
diagrams shown in figure 4, which for t ^ (rioA)^^ gives A = 2nQ to lowest order 
in a small density expansion [[7|. Aside from this term our action is now linear in 
the response fields, which we integrate out to yield the equations of motion: 

f + v,y^ = BV^^ - KM'^' - (22) 
^+„„,|^.flV^.^, (23) 

where Ae// is a new effective reaction rate constant, found by summing the diagrams 
shown in figure 5, giving 



in the limit v^t ^ 1. Note that this result ensures that the density amplitudes are 
Vo dependent even above the upper critical dimension dc = 2. If we now average 
equation ( ^2]) over the initial conditions, then we have 

since V^(0) = and dx{<f)) = 0. However we can see from a diagrammatic expansion 
for {(f)) (figure 7) that the only diagram contributing to the value of (ip'^) in (|25|) is 
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the single loop, which we now evaluate: 
(V^^) = ^ j dxdyd'^-^zG{x,y,z,tf 

dxdyd^^'^z ( 3 



A 



exp 



6(x - \vQytf 



(26) 



{4:'KDtY \{votY + 12 j \ Dt[{votf + 12] 2Dt 2Dt \ 



For t>ot ^ 1 and d = 2, this gives the result: 



(27) 



(2^ 



We can now find an upper bound solution to (|25|) by replacing (0^) by (0)^ (see 
[0 for a proof). Calling this upper bound solution /, we have for vot ^ 1 



df 



XeffAVS 



(29) 



dt ' IQnDvot^' 

This can be solved by making the substitution / = ii/Xu, giving the upper bound 

1 



(0)</ 



2A 



+ 



1 



\(2A 



eff 



AV3 
)2 ^ IGnDvo 



(30) 



m d = 2. However, we can also find a lower bound for (0) by noting that 0(x, t) > 
|'?/'(x, t)|, or equivalently that a(x, t) and 6(x, t) are everywhere non-negative. We 
can prove this rigorously using the effective field theory equations for a and b: 



da 2 d^ ^ 

- = DV a-v,y--Kjfah 



(31) 

effUU. (32) 

We now assume that the (smooth) fields a and h are initially everywhere non- 



96 9b ^ , 

- = 6 - voy^ - Kfjah. 



negative. Suppose, at a later time, a = at a point, then a > locally around the 
point, implying that it is a local minimum. Hence V^a > and dxa = 0, meaning 
that dta > 0. For a region of a = 0, then we have dta = inside the region. On 
its boundaries we have V^a > and d^a = 0, giving dta > at these points. As a 
result the fields cannot pass through zero and will remain non- negative. 
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Since we have (j) > {ipl, it follows that (0) > In addition, at long enough 

times, we expect to have a normal distribution - a result of ip satisfying a (mod- 
ified) diffusion equation. Consequently (as in |^), we have 

P[^(t)]ocexp|--^^|, (33) 



and therefore 



{<p) > mm = Jzmy) 



TT 



\ 



^^^(«)-', (34) 



for d = 2. However since {g'^) > (g)^ for any real g, then we also have 

(0 - < ((0 - 1^1)') = (0^) + (^') -2m\), (35) 

and using > IV'I this gives us (in c? = 2): 

(0 - < (0^) - (^') = - ~^(0) ~ Oit-'). (36) 

In other words, for d = 2, the bound on the corrections is of the same order as the 
density. Consequently we cannot say that the density asymptotically approaches 
the lower bound (as could be said for {no\)~^ <^ t <^ Vq^ and 2 < d < A), only 



that it lies somewhere between the limits supplied by ( pOD and (|3^. In this respect 
the case of c? = 2 and t ^ t ^ (noA)^^ is similar to that of c? = 4 and 
(noA)^^ <^ t <^ Vq^. In addition both situations retain the mean field decay 
rate, but with modified amplitudes, which for d = 2 and t ^ Vq^, t ^ (ngA)^^ 
depend on vq. 

Note that in the limit XIjj-A/vqD'^ ^ 1, the upper bound reduces to 

/ 



i.e. the upper and lower limits differ simply by a numerical factor of (2/7r) 2 . On the 
other hand, if AgjjA/-uo-D^ ^ 1, then / ^ (Ae//t)^^. In the limit of strong shear, 
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where any reaction zones are broken up almost immediately, we have recovered a 
mean field decay from this upper bound. 

Finally, we can see that the lower bound ( p4D is of limited usefulness in the small 
A or large vq limits, as the bound decreases with either increasing vq or decreasing 
A. However, we can use the fact that XIj^jA/vqD'^ is dimensionless in d = 2 
to obtain an improved expression for the densities, by performing a perturbation 
expansion with this parameter. From (|25|) it follows that the zeroth order term of 



this series is a constant, equal to the small Agr,A/t>o-D^ limit of the upper bound: 



(a) 



1 VSXjffA 



^efft 







2\ 




V 




) 





(3^ 



5 Numerical Results 

In order to confirm some of our theoretical predictions we have performed Monte- 
Carlo simulations in ci = 2. Initially a square lattice of size x Ly was populated 
with equal numbers of randomly distributed A and B particles. The evolution of 
this initial configuration was simulated using the rare event dynamics (RED) tech- 
nique (see also |@). In this Monte Carlo method, the time increment is determined 
by the current configuration: if many changes in the configuration are likely, then 
the time increment is small, whereas if the configuration is very stable, the time 
increment is large. In RED, a list is made of all possible changes to the configu- 
ration (events) together with the expected time after which each event will occur 
(rates). In the present model, two distinct types of events could occur: 

(1) A and B particles at (x, y) could hop to neighbouring sites with rates: 
UP: Rate= 1 

DOWN: Rate= 1 

LEFT: Rate= 1 + s 



RIGHT: Rate= 1 - s 



y—0.5Ly 
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where Vq = 2s/ Ly. Note that the possible values of the shear gradient Vq were 
restricted to ensure that the hop rates remained everywhere positive (i.e. s < 2). 

(2) Each A particle could react with each B particle on the same lattice site, 
with a reaction rate A. The simulations also allowed multiple occupation of each 
lattice site, in accordance with our theoretical description of the system. 

One step in a RED simulation consists of incrementing the time scale with 
At = 1/ X)j(rj), where Vj is the rate for event j, and then allowing selection (and 
execution) of an event. The probability that event i is chosen is equal to pi — 

For an efficient implementation the events are organized in a binary tree, where 
each branch contains one event and has a weight equal to the rate of that event. 
The weight of a parent node is equal to the sum of the weights of its children. 
As the root node contains the sum of all rates, the time increment At is easily 
obtained: it is the inverse of the weight of the root node. To select a particular 
event i with a probability proportional to its rate rj, we start in the root node, 
descend to either its left or its right child with a probability proportional to their 
weights, and iterate this process until we have reached the bottom of the tree. The 
selected event is then executed and the weights in the tree of all events whose rates 
have changed, plus their parents, grandparents, etc. are updated. The use of the 
binary tree assures that the CPU time required for one step in the RED simulation 
scales with the logarithm of the size of the tree. 

In our simulations we have adopted periodic boundary conditions in the direc- 
tion of the shear flow (in the x direction), but with hard wall boundary conditions 
perpendicular to it {at y — and y — Ly). 

The results of simulations on 800 x 200 and 900 x 300 lattices are shown in 
figures 8 and 9 respectively. Data for the shearless cases (s = 0) are averaged over 
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20 runs, whereas those with shear (s = 1) are averaged over 40 runs. At very 
early times the decay is fairly slow (especially in the no = 0.04 case) as it takes 
time for the particles to "find" each other. However at slightly later times we see 
a straight line for the densities, which is quite well described by the theoretically 
predicted decay. Notice that the simulations both with and without shear are 
indistinguishable at these early times, indicating that the shear really does have a 
negligible impact early on (as was predicted in section 4). However for s 7^ 0, we 
see a crossover at times of order Vq^ to a new regime, which is fairly well described 
by a decay. 

Notice that at late times the densities in the shearless simulations also tail off, 
and appear to fall away much more quickly than the predicted t^^/^ decay. Similar 
effects were also seen in the original paper by Toussaint and Wilczek [|I|], where 
the decay (for d < A) was first proposed. We interpret this result as being 
an exponentially fast density decay caused by finite size effects, which become 
appreciable when the segregated domain sizes become of the same order as the 
system size. At large enough times these segregated zones will have grown to a size 
where only two such regions exist (one A rich and the other B rich). The depletion 
zones of these segregated areas will typically be of the same order as the system 
length. Hence, if we have just two zones remaining in a system of size L x L, then 
the densities would decay according to: 

^(L'^a) ~ -JL ~ -(aD/L)L ~ -Da, (39) 
at 

where J is the flux of particles flowing into the reaction zone. This implies 

a ~ exp (- (const. )Dt/L2). (40) 

In order to test this prediction, we simulated an initially completely segregated 
system of size 200 x 200. The t = density profiles for each of the two species were 
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set up to rise linearly from zero (around y = Ly/2) to a maximum on opposite 
boundaries (at ?/ = and y = Ly). This configuration mimics the effects of 
having a diffusion length of the same order as the system size. The results of these 
simulations (averaged over 5 runs) are shown in figure 10, where an exponential 
decay is found, in agreement with the above theory. 

A further feature of the simulations concerns the effects of the hard wall bound- 
ary conditions. Segregated zones neighbouring the hard wall borders at the top and 
bottom of the system are likely to be more stable, as for these zones fewer reactions 
are occurring on their boundaries relative to similar regions in the bulk. Hence, if 
the system is roughly square, then we might expect the reaction zone at late times 
(when only two segregated regions remain) to lie roughly along y = Ly/2. This 
effect has indeed been seen in our simulations in the case of zero shear. 

Finally, in figures 11-13, we show the results of various "snapshots" of an 800 x 
200 system as it evolves in the presence of a shear fiow. The effects of the shear in 
tearing apart the reaction zones are clearly illustrated. 

6 Conclusion 

In this paper we have studied theoretically and numerically the effects of a linear 
shear fiow on the two species reaction A + B ^ ^. For t <^ Vq'^ we found that 
the system was essentially unaffected by the velocity fiow, and the species again 
segregated into A and B rich zones for d < A. In the regime (noA)^^ <^ t <^ Vq^ and 
with Dt ^ Uq'^^'^, our results for the densities (to leading order in a small density 
expansion) reproduced those of the shearless case 0. However at large times 
t ^ Vq^, the critical dimension for the onset of anomalous kinetics is reduced from 
dc = 4: to dc = 2. Consequently, for c? > 2, a mean field decay is adopted, though 
with a modified vq dependent amplitude. Notice, however, that our calculations are 
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for the case of equal diffusion constants. Although it would be possible to treat the 
case of Da 7^ Db, we believe (as was found in that such a modification would 
be qualitatively unimportant, provided both diffusion constants remain non-zero. 

To go beyond the calculations presented in this paper, one could attempt to 
generalise our analysis to the case of a non-linear shear fiow. Unfortunately, ex- 
cepting the very simple arguments of section 2, it is not clear how to make analytic 
progress in this the equation for the Green functions ( |12|) becomes very 

much harder to solve. However, it should prove possible to incorporate into our 
model a repulsive force between like particles. This would make our system more 



similar to those considered in ITT], [T^, |13|, where a same species hard core exclu- 
sion rule was imposed. In this way we could use the field theoretic formalism to 
investigate the effects of exclusion on reaction-diffusion systems with shear or drift. 
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Figure 1: Vertices given by the field theoretic action. 




Figure 2: The classical density {(f)) a, given by (a) the sum over tree diagrams, or 
(b) an integral equation. 




Figure 3: The diagrammatic sum for the response function. 

>o ■ X • 

Figure 4: The initial term A is generated by the sum of diagrams shown here. Only 
the leading order diagram is evaluated. 
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Figure 5: Diagrammatic expansion for the effective coupling constant. 



l(k,t) 




Figure 6: The sum of diagrams contributing to the primitively divergent vertex 
function X{k, t). 
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Figure 7: Diagrammatic expansion for (0), using the initial terms Uq and A. 
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Figure 8: Simulation results on an 800 x 200 system, with A = 1000, Uq — 0.1, and 
s = 1 (diamonds) or s = (plus signs). The straight lines have gradients of —0.5 
and —1 respectively. 
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Figure 9: Simulation results on a 900 x 300 system, with A = 1000, Uq — 0.04, and 
s = 1 (diamonds) or s = (plus signs). The straight lines have gradients of —0.5 
and —1 respectively. 
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Figure 10: Simulation results on a 200 x 200 system with complete initial segrega- 
tion and A = 1000, s = 0. 
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Figure 11: Snapshot of an 800 x 200 system, with s = 1, no = 0.1, A = 1000 at 
t = 300. 



□ □ ■ 1^ f H ° □ □ ° ■ ■ ' ■ ■ 

Figure 12: Snapshot of the same 800 x 200 system at t ^ 900. 




Figure 13: Snapshot of the same 800 x 200 system at t = 1500. 
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